pss header will be provided by the publisher 



Bottom-up approach to high-temperature superconductivity 

Hyekyung Won\ Stephan Haas^, and Kazumi Maki^ 



' Department of Physics, Hallym University, Chuncheon 200-702, South Korea 
Depa 
USA 



^ Department of Physics and Astronomy, University of Southern California, Los Angeles, CA 90089-0484 



Key words cuprate, d-wave superconductivity, d-wave density wave 
PACS 74.70.-b 

Since the discovery of high-temperature superconductivity in the cuprates a theoretical understanding of 
their phase diagram has remained one of the major outstanding problems in condensed matter physics. Here 
we propose an effective low-energy Hamiltonian which produces both d-wave density wave (dDW) and 
d-wave superconducting (dSC) solutions within the BCS mean-field theory. This model predicts that (a) 
the observed pseudogap phase is a dDW state, (b) the superconducting phase is a d-wave BCS state, and 
(c) in the underdoped regime there is a gossamer superconducting state, i.e. dSC in coexistence with dDW. 
Moreover, this theory naturally explains the Uemura relation, the reduction of the quasiparticle density of 
states at the Fermi level, and the salient features in the tunneling conductivity measured in underdoped 
Bi2212. 
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1 Introduction 

In 1986, the discovery of high-temperature superconductivity in the cuprates took the physics community 
by surprise. Ql This was the starting point of a new era in condensed matter physics. Many practical appli- 
cations of these new compounds were envisioned, and at the same time an intense debate arose regarding 
the origin and possible mechanisms leading to this new phenomenon. Enz recorded the often confusing 
discussions of the early days of high-Tc research in his beautiful textbook. One of the most influential 
contributions to the theory of these materials was provided by Anderson's "dogmas". f3l He stated that the 
cuprate high-Tc phase diagram arises from an inherent competition between a Mott insulator phase and 
s-wave BCS superconductivity in these materials. In order to model high-Tc superconductivity, he pro- 
posed a two-dimensional one-band Hubbard model in combination with a resonant valence bond (RVB) 
wave function. A great portion of the theoretical community in the field has since embraced these dog- 
mas. However, unfortunately we still remain without a clear vision as to where these dogmas are leading 
us.|4| Around 1990, Scalapino and othersfsl pointed out that a perturbative analysis of the 2D Hubbard 
model in the weak-coupling limit produces d-wave superconductivity. Indeed, a d-wave superconducting 
order parameter was experimentally established around 1994 for single crystal samples of optimally doped 
Bi2212, YBCO and LSCO, using powerful angle resolved photoemission spectroscopy (ARPES)f6l and 
elegant Josephson interferometry[X^. These observations motivated us to investigate d-wave supercon- 
ductivity within the BCS framework. l9l [Tol [TTl 

Before elaborating further, let us first examine the generic phase diagram of the hole doped high-Tc 
cuprate superconductors shown in Fig. 1. From the beginning, this phase diagram has been hotly debated. 
Around the year 2000, a few groups suggested that the pseudogap region can be described by a d-wave 
density wave (dDW) phase. Indeed, the giant Nernst effect observed in the underdoped Bi2212, YBCO and 
LSCOl 17, 18 19 1 and the angle dependent magnetoresistance in Yo.68Pro.32Cu04 |20| have been found 
to be fully consistent with dDW. 12 II 1221 In past work, we have shown that these are consequences of the 
Landau quantization of the quasiparticle spectrum in a magnetic field, analogous to earlier considerations 
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Fig. 1 The phase diagram for the high-Tc cuprates. p denotes the hole doping concentration. PG is the pseudogap 
region. 

by Nersesyan et al. 1231 1241 Moreover, we note that the Fermi arcs (or pockets) in the (tt, tt) directions, 
observed by ARPES, follow directly from dDW. 1161 l25t 26 1 Furthermore, it is by now well established 
that the overdoped regions of the cuprates can also be described in terms of a d-wave BCS model. LI llll5l 

Recently, Laughlin|27| pointed out that the Gutzwiller operator which is commonly used in the RVB 
wave function is not mathematically tractable, and proposed to replace it by a less constrained Jastrow 
operator. He named the resulting coexistence phase "gossamer superconductivity", i.e. a condensate with 
a reduced superfluid density and reduced density of states near the Fermi surface. Such a reduction of the 
quasiparticle density of states has recently been observed by Tallon et al l28ll29l by means of a thermody- 
namic analysis and the effect of Zn impurities in YBCO over a wide doping range. In section 4 we will 
return to the characterization of gossamer superconductivity. Note also that gossamer superconductivity 
emerges naturally from the phase diagram in Fig. 1 as a coexistent phase of d-wave superconductivity 
in the presence of a d-wave density wave. 1 141 1151 1161 In the following we will explore this gossamer 
superconductivity phenomenon in detail. 



2 Effective Hamiltonian 

In this section, we construct an effective low-energy Hamiltonian that constitutes the basis of the bottom- 
up approach. This approach should be viewed as an alternative to the common top-down approaches 
originating from higher-energy Hamiltonians, such as the t-J and Hubbard models. Considering the energy 
scales of these models, e.g. typically a Hubbard U of the order of lO^K, it turns out to be a rather difficult 
task to arrive at superconducting phenomena that exist at scales of Tc'^ lO'^K. Based on the renormalization 
group analysis of 2D electron systems|30| we understand that the normal state is a Fermi liquid, i.e. not 
a Luttinger liquid or bosonic liquid. Here we define the Fermi liquid via a quasiparticle Green function 
which has simple poles, a definition that is consistent with Shankarf301 and Landau|31 1. Furthermore, we 
know that at sufficiently low temperatures the normal state becomes unstable against infrared divergences 
in the 2-particle and/or 2-hole channels, implying superconductivity, or unstable against divergences in the 
particle-hole channel, implying density wave phases. 
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The effective low-energy Hamiltonian for such a system is given bv ll4l 1321 1331 

(A2(fc)<^cL,^^ + A;(fc)c_fe,xCfe,T) - .9r'|Ai(fc)r - 52"'|A2(fc)|^ (D 



fc 



where the ampUtude and angular parts of the order parameters separate via Ai(fc) = Ai/(fc) and A2(fc) = 
A2/(fc), and Q ^ (tt, tt) is the nesting vector Two self-consistent gap equations follow directly from this 
Hamiltonian: 

= ^7%)yE^''^/W(4+Q..c,,.)> (2) 

Here Ai and A2 are the order parameters of dDW and dSC respectively. In the following, we use 
f{k) = cos(2(/)) as the angular dependence. A similar Hamiltonian has been considered in related work by 
Thalmeier [32] However this study was limited to conventional DW and conventional SC, and to the case 
of vanishing chemical potential /i. As we shall see here, /i is an important control parameter in the present 
model. fTBil Moreover, when both DW and SC are conventional, there is little room for their coexistence; in- 
stead phase separation is the rule. 1 32 1 On the other hand, when both order parameters are unconventional, 
there is ample opportunity for their coexistence. L33J This fact will be exploited in the following. 
The Nambu-Gorkov Green function l34l of this model is given by 

G'~^(fc,u;„) = iun - efcp3(T3 + fias + | Ai | exp(-i0ip3)/(fc)pi(T3 + IA2I exp(-i02O-3)/(fc)o-i, (4) 

and the corresponding spinor field is 

= (cl^^,C^k,-a, cl+Q,a^C^k-Q,-a) ■ (5) 

One notes that unlike in Ref. 1141 . the present G^^{k, ojn) possesses two Abelian gauge transformations 
associated with cpi (sliding motion of dDW) and (j)2 (supercurrent in dSC). The determinant of G^^{k, Un) 
is given by 

D = det\G-\k,LU,,)\ - {ujI + el + fi^ + |Ai(fc)p + |A2(fc)|2)' - 4^^ (,2 ^ |Ai(fc)|2) . (6) 
Using this result, the quasiparticle energy is found to be 




E^±^^^el + \A^{k)\^T^^) +|A2(A:)|2, (7) 
which agrees with earlier results. 1121 1301 Finally, the quasiparticle density of states is given by 



N{E) 
No 



E a 

Re , IT ^ 



^ei-Aif{kYTlA -Aifik) 



^E^ - A'ifiky 



(8) 



where + and - stand for the positive-energy and negative-energy solutions respectively. 

In Fig. 2 the quasiparticle density of states is shown for a particular set of parameters. For these 
parameters, we observe clear dips in the vicinity of E=0 , as well as a quasi-linear dependence on energy. 
In the regime Ai — /i < A2 one of the peaks splits into two peaks. Although such a split peak has not yet 
been observed experimentally, except in the presence of Ni impurities, 1 35 1 the present result is consistent 
with recent measurements on underdoped Bi2212. l36l Moreoever, this result is rather different from earlier 
work by Zeyher and Greco ll37l who assumed /i = 0. 
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Fig. 2 Quasiparticle density of states in a gossamer superconductor with the energy scale set by the superconducting 
energy gap A2 = 1. 



3 D-Wave Density Wave Phase 

Equations 2 and 3 can be transformed to|i33l 

Ar' = 47rT^i?e(/(fc)2d-i) 

n 



(9) 



(10) 



where d = \{\/ + A|/(fc)2 - ^lf + A^/(fc)2] V2, Here Ai = .giTVo and A2 52A'^o are dimensionless 
coupling constants, and A^'o — ./V(O) is the quasiparticle density of states in the normal state. 

First, the phase diagram of the pure dDW states is easily obtained by setting A2 — 0. I1I6J In this case, 
Eq. 10 reduces in the limit Ai — *■ to 



(11) 



where Td is the transition temperature for dDW and ip{z) is the di-gamma function. Using Tcio=800K 
1381 . one arrives at the phase diagram shown in Fig. 3. Note that Eq. 13 is the same for s-wave and d-wave 
superconductors 1 39 , 40 , 4 1 1 in the limit when the Pauli term dominates over the orbital term. 

It is observed that Td bends backwards in the region /i/Aio < 0.558. A similar diagram has also been 
found in Ref. \37\. However, if we additionally allow a spatial variation of the dDW order parameter Ai, 
we obtain 



In 



i?e ( (1 ± cos(40)) 1/; ( 



(12) 



where p — v\q\/2ij.. This yields the extended portion shown in Fig. 3, analogous to the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) 1421 1431 state in d-wave superconductors.! 44] There is a further transition 
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Fig. 3 The phase diagram of the high-Tc cuprates. The shaded area denotes Tc^il^), whose determination from the 
present model will be presented in a future publication. 

when the q-vector is rotated from the [100] to the [110] direction. Finally, the dDW regime terminates 
when /i/Aio = 0.824. We call these periodic dDW phases dDWII and dDWIII respectively. 



4 Gossamer Superconductivity 

Now we can ask how dSC appears on top of this dDW background. In the following, we shall limit 
ourselves to the region Ai A2, /i. Then we can deduce in the vicinity of the superconducting transition 
temperature, Tc2, the quasiparticle density of states at the Fermi surface which is given by 



g{fi,0) = (g(0,fc)) = -xK{x), 



(13) 



with X = ^/Ai{fi) and K{z) is the complete elliptic function of the first kind. g{fJ.,0) resembles the 
quasiparticle density of states deduced from the analysis of Zn-impurities L28J . The corresponding low- 
temperature entropy is given bv ll29l 



T 



27C(3) T 



27r2 



A2(m)' 



(14) 



where = tt'^Nq/3. On the other hand, the superconducting transition temperature and free energy are 
controlled by 



as 



4x 

giifi, 0) = 2(cos2(2</))g(0, k)) = —{K{x) - E{x)), 

TT 



T.ifi) = 1.136Ai(Ai)exp[-(Ar^-A2-i)5-i(A*,0)], 
Uo = -\Nogi{fi)[A2{fi,0)f. 



(15) 

(16) 
(17) 
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Eq. (16) appears to somewhat overestimate Tc(/i), and hence a more detailed treatment needs to be devel- 
oped in order to be more realistic. The functions Ai (/i) / Aiq, g{iJ.,0) and gi {fi, 0) are shown in Fig. 4. The 
dependence of the functions g(/i,0) and gi{fi,0) on /i/Aio is in agreement with available experimental 
data. l28ll29l 




Fig. 4 Functions Ai (/i) / Aio, g{jJ., 0) and gi {/j,, 0), controlling the low-doping region of the high-Tc superconduc- 
tors. 



By solving the coupled gap equations in the regime A2, /i ^ Ai we find ll6l 

P'^^> - - A2(0) ' ^^^^ 
_ 1 AiO.,0) 

" 21n2 A2(0) ' ^ ^ 

A-20,,0) - —rPPsio). (20) 

TO* 

These are essentially the Uemura relations. [45] Therefore, if we limit ourselves to the deeply underdoped 
region, many experimentally observed features of gossamer superconductivity follow naturally from the 
present model. 



5 Concluding Remarks 

The model treated here is based on an effective low-energy Hamiltonian which describes dDW and dSC 
states with the chemical potential as a control parameter This theory accounts for the following principle 
features of high-temperature cuprate superconductors: (a) the normal state is a Fermi liquid, (b) the pseu- 
dogap phase is a dDW (more recently a d-wave spin density wave has also been suggested ll46ll47l '). (c) the 
superconductivity in the optimal to overdoped regime has a BCS d-wave order parameter, and (d) in the 
underdoped regime there is gossamer superconductivity, i.e. dSC coexisting with dDW. 

Recent related studies on heavy fermion compounds, such as CeCoIn5 under pressure] 48j, and organic 
conductors, such as /3"-(BEDT-TTF)4[N30)M(C204]C5H5N with M = Ga and Cr k49ji50i, have revealed 
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many parallels between the high-Tc cuprates and these systems. These include (a) a layered structure or 
quasi-two-dimensionality, (b) d-wave superconductivityf5T32|, and (c) d-wave density wave phases |47| 
|53, J4 1. In the heavy fermion and organic conductors the horizontal axis in Fig. 3 needs to be replaced by 
the external pressure P, but otherwise their phase diagrams look very similar. This suggest strongly that the 
present model is rather universal and applies to many strongly correlated electron systems. 

Let us finally note that gossamer superconductivity is not necessarily restricted to dSC and dDW. For 
example, recent experiments on the Bechgaard salts (TMTSF)2PF6 at ambient pressure lSSl 1561 1571 and 
measurements of the angle dependent magnetothermal conductivity in URu2Si2 1 58 1 suggest that there are 
other kinds of gossamer superconductivity, e.g. f-wave superconductivity coexisting with a d-wave spin 
density wave. Hence it will not be surprising if similar coexistence states will soon be discovered in related 
strongly correlated electron systems. 
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